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1 Introduction 


The measurement of the Higgs boson mass at the Large Hadron Collider [Tj not only fully 
confirms the validity of the Standard Model (SM) around the electroweak scale, but also 
opens a possibility for a precise study of the applicability of the SM at energies of the 
order of the Planck mass. Renormalization group (RG) equations, which determine the 
dependence on the renormalization scale /i of the running parameters, which are usually 
defined in the modified minimal-subtraction (MS) scheme of dimensional regularization, 
play an essential role in such analyses. In the SM, the corresponding RG functions, whose 
knowledge had been limited to the two-loop order for a very long time [3], have recently 
been evaluated at three loops for all running parameters, including gauge, Yukawa, and 
scalar self couplings [3j. 

The other aspect of the problem is the matching between the MS parameters and the 
physical observables, which gives rise to threshold correct ions @ These not only include 
terms of the form In /i 1 2 , but also non-logarithmic ones. The relationships between the 
MS and pole masses of the intermediate bosons were obtained at the two-loop level in 
Refs. [H [5]. As for the threshold corrections to the top and bottom quark masses and 
Yukawa couplings, the situation is as follows. The QCD corrections, which are dominant, 
are available at one [6], two m, three [9], and four [10] loops. The two-loop result in the 
supersymmetric extension of QCD was obtained in Ref. HU. The one-loop electroweak 
corrections, of order 0(a), were first considered in Ref. [13] . The two-loop mixed 0(aa s ) 
corrections were provided for the bottom quark in Ref. [13] and for the top quark in 
Ref. [HJ[15][l6]. Recently, also the two-loop electroweak corrections of order 0(a 2 ) have 
been obtained in the gaugeless-limit approximation im Also the threshold corrections to 
the self-coupling constant of the scalar field were intensively studied in the literature. The 
0(a) corrections were evaluated a long time ago in Ref. pT8] and the 0(aa s ) ones recently 
in Ref. [16] . As for the 0(a 2 ) corrections, the leading term was found in Ref. [19], and an 
interpolation formula, which also includes subleading contributions, was given in Ref. [20], 
These analyses were recently revisited in Ref. pT] by providing precise numerical results. 

In this paper, we systematically present the complete two-loop threshold corrections, 
from the orders 0(a), 0(aa s ), and 0(a 2 ), to all the running parameters of the SM 
independently obtained by an analytic calculation. This includes the masses of the W, Z, 
and Higgs bosons (mw, mz, i^h) and those of the top and bottom quarks ( m t , m&) as well 
as the gauge couplings ( g , g'), the Higgs self-coupling (A), and the top and bottom Yukawa 
couplings ( y t ,yb )• In contrast to Refs. [20] [2TJ E22] > all our calculations are performed in 
R^ gauge keeping the gauge-fixing parameters free, which allows us to explicitly track the 
£ dependencies and so to ensure the gauge independence of the threshold corrections and 
the MS parameters. The tadpole diagrams turn out to play a crucial role in this (see 
Subsection 12.21) . 

This paper is organized as follows. In Section [2] we set up the stage for our calculation 
of the threshold corrections. In Subsections I2.1H2.41 we discuss the various ingredients 
entering our analysis. Our results are presented in Section [3] and Appendix [A] In Ap¬ 
pendix [B] we also list the MS renormalization constant of the Higgs boson mass. 

1 The usage of the term threshold corrections in this context is to indicate that the initial conditions for 

the RG evolution of the MS parameters are determined at some low-lying scale. This term also appears 
in the effective-field-theory language, where it carries a different meaning. 
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2 Setup 


The SM may exist in two different phases: the symmetric phase and the phase with the 
spontaneously broken symmetry. The phase is determined by the potential of the scalar 
held </>, 

V{4>) AO f 0) 2 , (1) 

where is a mass parameter and A is the self-coupling constant of the scalar field. 
While stability requires A > 0, the term m 2 ^ can be either positive (symmetric phase) or 
negative (broken phase). In the symmetric phase, the SM is naturally parametrized by 
the following set of parameters: 

g,g',X,m^,y f . ( 2 ) 

In the broken phase, it is convenient to choose an alternative set of parameters, namely 


e,m w ,m z ,m H ,mf, 


(3) 


where e is the electromagnetic gauge coupling. 
At the tree level, the parameters in Eqs. (J2| 


and © can be related to each other by 


g 9 


1 

I2 


(4) 


for the couplings and 


4 rrtfy 


= 9 


4m| 


= g 2 + g ' 2 


m 


H 


2v 2 


= A, 


2mj 


y) 


(5) 


for the masses, where we have introduced the vacuum expectation value of the scalar field 
v ~ 246 GeV, which characterizes the broken phase. This parameter is not independent, 
since the parameter sets in Eqs. (J2J) and © each fully determine the theoryH In fact, we 
have 


-m. 


v = 


in terms of Eq. (ED and 


( 6 ) 


(7) 


v 2 4m^ (1 - m^/m|) 

in terms of Eq. (ED . 

Equations (ED-Q are subject to radiative corrections. However, we have the freedom 
to choose 16 out of the 32 parameters in Eqs. (ED and (ED to be “observables” and to con¬ 
sider the remaining ones to be derived parameters of the theory. At energies of the order 
of v, the parameters in Eq. (ED are intuitively closer to what we would call “observables,” 
especially if we define them in the on-shell renormalization scheme. Then, these include 
Sommerfeld’s fine-structure constant ctTh = e^ h /(47 t) as measured in Thomson scattering 
and the pole masses M\y, M z , Mh, and Mf, which are the zeroes of the inverses of the 
propagators of the respective particles. By contrast, the parameters in Eq. (ED are more 
suitable for studies of the RG evolution. In this work, we shall discuss the relationships 
between these parameters and the radiative corrections to Eqs. 0 - 0 . 


throughout this paper, we take the quark mixing matrix to be unity. 
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As usual, we consider Eqs. (HJ—(JTJ) to be valid through all orders of perturbation theory 
by definition. For the reason explained above, we define the parameters in Eq. (J2]l, which 
appear on the right-hand sides of Eqs. © and ([5]), in the MS renormalization scheme, 
which then carries over to the parameters in Eq. ([3]), which appear on the left-hand sides 
of Eqs. (HD and © (see Subsections 12. 3l and l2.4p . The /i dependencies of the parameters in 
Eq. © are governed by the RG equations. These allow us to run the parameters from a 
few GeV way up to the Planck mass. The relevant (3 functions are available at two j2] and 
three loops |3]. The masses in Eq. © also obey RG equations. The corresponding RG 
functions are available at the two-loop order from Refs. mum and are partly recovered 
here by an independent calculation. As explained in Refs. mu, we can always relate the 
RG functions of the unbroken and broken phases, which serves as a welcome check for the 
correctness of our results. 

As for the initial conditions for the RG evolution, the MS parameters are usually 
expressed in terms of the parameters of the on-shell renormalization scheme, «Th ; ^w, 
Mz, Mh , and Mf. The matching scale where this is done is typically chosen to be of 
the order of Mz or M t . In this case, one may safely neglect the masses of all fermions, 
except for the one of the top quark and possibly also the one the bottom quark. We shall 
retain the full dependence on the bottom-quark mass at one loop, but neglect it at two 
loops. The strong-coupling constant a s (n) is always defined in the MS renormalization 
scheme. Finally, v is related to e via Eq. ( 0 . Phenomenologically, it is more convenient 
to express the MS parameter v{/x) in terms of the Fermi constant Gf, which is measured 
in low-energy processes of the weak interaction, such as muon decay, through the exact 
relationship 


1 + A r(n) 
u 2 (/i) 


2 1 ' 2 G F = (8) 

(see Subsection im 


where Ar(/i) is an appropriate variant of Sirlin’s Ar parameter |) 

Inserting Eq. © in Eq. © and accommodating the relationships between the MS and 
pole masses, one may cast the relations of interest in the form 


g 2 (g) 

= 2^ 2 Gj?M^[1 + 5w(g)\ , 

(9) 

g 2 (g) + g ,2 (g) 

= 2 5 ^ 2 GpM^[l + Sz(g)\ , 

(10) 

A(/i) 

= 2- 1 / 2 G f M 2 h [1+Sh(v)}, 

(11) 

vM 

= 2^G l J 2 M f [l + 5 f (v)]. 

(12) 


In turn, substituting Eq. © in the right-hand sides of Eqs. ©- (1I2|) and equating the out¬ 
come with the right-hand sides of Eq. ©, we may write the MS to pole mass relationships 
as 


ml(fi) = M^[l + Af(fx)][l + S x (n)], x = W : Z,H, 

m x {n ) = M x [l + Ar(/i)] 1/,2 [l + <5 X (//)], x = t,b , (13) 

which are valid to all orders of perturbation theory. It is the goal of this work to evaluate 
Ar(/i) and 5 x (n) for x = W, Z, H, t, b through order 0(a 2 ). 


2.1 Vacuum expectation value and Ar(p,) 

We start our discussion with the evaluation of the threshold corrections to u(/i), i.e., of 
Ar(/i) in Eq. ©. As already mentioned above, this quantity is very similar to Ar in 
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the on-shell scheme [25] . which contains the radiative corrections to the muon lifetime 
that the SM introduces on top of the electromagnetic corrections in the Fermi model of 
four-fermion interactions. The difference between Ar and A r{/i) is that the former is 
defined in the on-shell scheme according to 


G f = 


7r«Th 

V2M*y{ 1 - il^/Mf) 


(1 +Ar), 


(14) 


while the latter obeys the analogous relation with the on-shell parameters replaced by 
their MS counterparts, 


G f = 


7rct(/i) 

V2'm 2 w (n)[l - m^(/i)/m|(/i)] 


[1 + Ar(/i)]. 


(15) 


Consequently, the calculations of Ar and Ar(/i) are very similar. In particular, they are 
both based on the matching of the Fermi model and the SM and exhibit a factorization 
of the low-energy scales. The contribution of the light (massless) fermions to A r(/i) may 
be found in Ref. [24]. The MS parameter a(/i) deserves a detailed discussion, which will 
be presented in Subsection 12.41 

The framework for the evaluation of Ar or Ar(/i) at any loop order was established 
in Refs. [251 [25] • One may start from any amplitude involving the exchange of the weak 
charged current in the SM, e.g., A(e + v e —> n + u^), where e, u e , n, and are the electron, 
the electron neutrino, the muon, and the muon neutrino, respectively. As demonstrated in 
Ref. [25], there is a factorization theorem that allows for the convenient separation of the 
soft scales (see also the discussion in Ref. [27]). Applying this to the evaluation of A 
we only need to perform the renormalization procedure in the MS scheme. According to 
Ref. [25], we may write 


8m^(l - m%y/m? z ) 


(1 + Ar) — \ /^2,e^2,i/ e -^2,/i-^2,t/ M A(e + V e —> h + Vfi) 


hard 


(16) 


where Z 2 j is the wave function renormalization constant of the left-handed held com¬ 
ponent of fermion / in the MS scheme. The subscript hard in Eq. (Tl6l) implies that all 
external four-momenta and all the light-fermion masses are identically put to zero be¬ 
fore the integration over the loop momentaObviously, such a procedure generates a 
lot of infrared divergences in the calculation. These are regularized by the dimensional- 
regularization parameter e and cancel on the right-hand side of Eq. (ITFll . This cancellation 
is nontrivial, and this is exactly the statement of the factorization theorem [25] . Following 
this procedure, the evaluation of Ar is completely reduced to vacuum diagrams with one 
or two loops. In Appendix [A] we shall present an analytical expression for Ar in terms 
of MS couplings and pole masses. 


2.2 Role of tadpoles 

An important comment is in order here. In a theory with a broken gauge symmetry, it is 
important to take into account the tadpole diagrams. One- and two-loop tadpole diagrams 
contributing to the propagator of a particle are shown in Fig.Ul All such tadpole insertions 

3 In other words, all the loop momenta k^ are hard compared to the low-energy scales, such as the 
muon mass to m , that is \ki\ to m . 
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must be made not only in the Feynman diagrams contributing to the counterterms, but 
also in all the proper Feynman diagrams to ensure the gauge independence of renormalized 
scattering amplitudes. In particular, the mass counterterms are gauge dependent unless 
the tadpole contributions are included, as has been known for a long time [28]. This is 
also true for the threshold corrections, as was first observed for the one-loop electroweak 
threshold corrections to the Yukawa [T2] and Higgs self-couplings [IS] . 







Figure 1: One- and two-loop tadpole contribution to the propagator of a particle. H 
stands for the Higgs boson propagator with zero momentum transfer. 

In this work, we perform all calculations at two loops in R g gauge with four independent 
gauge parameters, related to the W and Z bosons, the photon, and the gluon. We verify 
by explicit calculation that the threshold corrections A r in Eq. (J8]) and 5 X in Eqs. (T9l) (TT2|) 
are gauge independent. This serves as strong check for the correctness of our results. 

Tadpole contributions are singular in the limit Mh —> 0. The most divergent terms 
scale as 1 /Mjj at one loop and as 1 /Mjj at two loops. This behavior is fully inherited by 
A r and the MS masses. However, nontrivial cancellations between A r and the MS masses 
take place in the threshold corrections &z, $t, and 5b, which render the latter finite 
in the limit Mh —> 0, i.e., also terms proportional to are canceled. As for the 

Yukawa couplings, this was noticed at 0(a) in Ref. [T2] and at 0(aa s ) in Ref. [13] for 
the bottom quark and in Ref. pi6] for the top quark. Here, we investigate the behavior 
for Mh —> 0 at 0(a 2 ) and find that Ar and the MS masses contain terms proportional 
to 1/Mfj, while 5w, 5z, 5 t , and 5b are finite in this limit. 

The situation is different in the Higgs sector. Already at one loop, the threshold 
corrections to m# and A contain terms that diverge as 1 /Mfj for Mh —> 0 [18j . At 
two loops, however, the 1 /Mjj terms cancel, so that the leading small- Mh behaviors go 
unchanged. This cancellation ensures that A stays finite in the limit Mh —> 0, since 
A ~ Mjj according to Eq. (TTT1) . 

The tadpole contributions may be quite sizable numerically. In Ref. mm , it was 
noted that one-loop electroweak correction to the MS to pole mass relationship of the 
top quark roughly compensates the QCD one. In particular, the tadpole contribution 
contains terms that are enhanced as N C M^/(M^Mjj), where N c = 3 is the number of 
quark colors, at one loop. At two loops, such terms appear in square. 

2.3 MS masses 

The pole mass M of a particle is defined to be the position of the pole of its propagator, 
i.e., the zero of its inverse propagator 0 For a scalar boson, such as the Higgs boson, one 

4 If the particle is unstable, then the pole position has a complex value. The latter is usually 
parametrized as p 2 ole = M 2 — iYM for bosons and as i> vo \ e = M — iT /2 for fermions, where M and 
T are the real pole mass and the total decay width of the particle, respectively. 
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thus needs to solve the equation 

0 = P 2 ~ m 2 H o - II hh(p 2 ) , (17) 

where p is the four-momentum, mu,o is the bare mass, and Hhh(p 2 ) is the self-energy 
function of the Higgs boson. The solution of Eq. (ITT)) gives the pole mass Mh as a function 
of rriH.o and other parameters of the SM. For the W boson, the appropriate equation is 
similar to Eq. mm, namely 


0 — p 2 ~ m w,o ~ bl ww,t{p 2 ) , (18) 

where n ww,t(p 2 ) is Hie transverse part of the TH-boson self-energy. The Z-boson case is 
somewhat more complicated because of the 'y-Z mixing. Diagonalizing the corresponding 
propagator matrix, one obtains 


o = p 2 - m| 0 - n zzAv 2 ) 


gj zjAf)_ 

P 2 - n 77 iT (p 2 ) ’ 


(19) 


where n zz,t(p 2 ), n 7 z,r(p 2 ), and U llT (p 2 ) are the respective transverse self-energy func¬ 
tions. Finally, for a fermion /, one has 


0 = ^ — m /)0 -£/($, (20) 

where Eis the self-energy of /0 The solutions of Eqs. (fT7]) - (l20l) may be found order 
by order in perturbation theory, by substituting the ansatze p 2 = m^l + Ki + + ...) 

and j> = mo(l + Hi + + • • ■), with k* being the i-loop corrections, in Eqs. (fT7l) - (ll9l) and 

Eq. (120)) . respectively. 

Alternatively, we may perform the mass renormalization in the MS scheme. The 
relation between the MS mass m(/i) and m 0 has the simple form 

m 2 0 = m 2 (p)(l + —+ — + ■■■) , (21) 


where the expansion parameter £ = (4 — d)/2 measures the deviation of the space-time 
dimension d from 4. The mass renormalization constants Z (j> have double expansions in 
the weak and strong gauge couplings, g(p) and g s {p), respectively. For the purposes of 
our two-loop analysis, we need to include the following terms 


Z {j) 



jLjLzU) , 

16vr 2 16vr 2 



Z^ + 


( 22 ) 


As already mentioned in Subsection 12.21 all the relevant mass renormalization constants 
are gauge independent upon the inclusion of the tadpole contributions. Explicit expres¬ 
sions through order 0(a 2 ) may be found for the W and Z bosons in Refs. mgs for the 
top and bottom quarks in Ref. ra, and for the Higgs boson in Appendix [Bl of this paper. 

5 Here, fermion mixing is neglected. The renormalization for mixed systems of spin-1/2 fermions was 
elaborated in Ref. [29]. 

6 The expressions for Zad 2 J of the W and Z bosons in Eq. (4.41) of Ref. [5] contain several misprints, 
which are corrected in Footnote 9 of Ref. [16]. 
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In this paper, we shall evaluate the MS masses m x (n) for x = W, Z, H,t,b through 
order 0(a 2 ). Besides gauge independence, also the fi dependence, which is dictated by 
the RG, provides a strong check for the correctness of our results. In fact, the full /i 
dependence of m x (n ) at two loops may be retrieved from the one-loop result for m x (n) 
and its two-loop anomalous dimension. To simplify the discussion in the remainder of 
this section, we shall only allow for one coupling constant, a(/i). The generalization to 
the case under consideration in this paper is straightforward. 

Adopting standard notations, the generic RG equations for a(/a) and m(/i) read 


da d lnm 2 

d In /i 2 ’ d In /a 2 


(23) 


The /3 and 7 functions in Eq. (123|) have perturbative expansions of the forms /3 = a 2 /3i + 
a 3 /?2 + • • • and 7 = 071 + a 2 72 + • • ■. In the MS scheme, R are just numbers, while 7 j 
are, in general, polynomials in squared-mass ratios. Expressing in turn the MS masses in 
terms of pole ones, we find the expression for 7 to assume the form 


7 — a(/r)R + cr(/i)(A 2 ln/x 2 + R) + ■ • • , 


(24) 


where R and A* are independent of /a. On the other hand, the MS to pole mass relationship 
has the form 


m 2 (fi) — M 2 [ 1 + a(/r)(Ai In /i 2 + B\) + a 2 (/i)(A 2 In 2 /i 2 + A 2 In /1 + B 2 ) + •••], (25) 

where A i: A', R, etc., do not depend on ji. The coefficients A 1} A 2 , and A! 2 in front of 
In"’ /i 2 with n = 1,2,... may be derived from the RG. Substituting Eqs. (1241) and (|25|) in 
Eq. (f23ll and comparing the coefficients of a m ln" /i 2 , we obtain 


Ai = R, 

2A 2 = T 2 + A 2 -RR, 

A' 2 = V 2 + BiTi — fiiBi. (26) 


Consequently, we may independently evaluate the coefficients of In n fi 2 with n — 1,2,... 
from one-loop results and the mass anomalous dimension 7 . In the MS scheme, the latter 
is completely determined by the single 1/e pole in the corresponding renormalization 
constant in Eq. (l!?T]h Taking into account both the weak and strong coupling constants, 
we have 


7 = 


2 dg 2 dg. j 


(27) 


On the other hand, 7 may also be related to the anomalous dimension 7 q f of the Fermi 
constant Gp and the RG functions in the unbroken phase of the SM |4j [5j 30] , 


2.4 MS fine-structure constant 


In this section, we discuss different definitions of Sommerfcld’s fine-structure constant 
a = e 2 /(4(/r) and their relationships to the Fermi constant GR- 

The relationship between the MS quantity cx(/i) and the on-shell quantity «Th defined 
in the Thompson limit is usually written as [3T] 


cx(/i) 


«Th 

1 — Ack(/i) ’ 


(28) 







where, to first approximation, Act(/i) is expressed through the vacuum polarization func¬ 
tion of the photon II'(0). The latter quantity is known to be subject to nonperturbative 
effects due to the hadronic content of the photon. The incorporation of the perturba¬ 
tive 0(aa s ) corrections was explained in Ref. |32j . The purely leptonic contributions are 
known through four loops [33]. A systematic discussion of higher-order corrections may 
be found in Ref. [31]. Through order 0(aa s ), we liav^l 


Aa(/i) = A a$ d (M z ) + 


Oi Th 
47T 


7 L 


w 


E 

Z=e,/i,r 


NcQ; 


a 


——L t + Cp- — {—L t + 15) 


47T 


+ n c Y, Q 


q=u,d,s,c,b 


--L z + — + C f — 

3 9 4tt 


-4 Li 


16Cs + 


55 


,(29) 


where A a^ d (M z ) = 0.02772(10) [35] is the hadronic contribution to II!)(0) and L x = 
In(Mj//r 2 ) with M x being pole masses. The four terms inside the curly brackets in 
Eq. (1291) are the bosonic, leptonic, top-quark, and perturbative light-quark contributions, 
respectively. To the same accuracy as in Eq. (j29]h we may also write the /i dependence 
of as 




a(M z ) 

a(M z ) f 11 80 a s \ M% 

47t \ 3 3 47t/ fi 2 


(30) 


Using a^(M z ) = 127.940(14) [35] as input in Eq. (j30j) . we obtain = 127.565 at 

order 0(a) and a _1 (M t ) = 127.540 at order 0(aa s ). The QED j3 function is available 
through Eve loops including QCD corrections [36] . 

Eqs. (1281) and (1291) relate a(M z ), which belongs to the core of electroweak physics, 
to the low-energy-QED parameter h, which is insensitive to the weak interaction due 
to the decoupling of the heavy particles. While this approach allows one to conveniently 
implement the RG and so to consistently accommodate all the QCD contributions of the 
orders 0(aa™ ), there are obvious drawbacks if one wishes to go beyond this approximation. 
First, it is less appropriate for the incorporation of electroweak physics, which comes into 
play by the exchange of W and Z bosons in the loops. Second, the definition of a(/x) 
from the QED two-point function alone has only restricted meaning. To relate a(/x) 
to electroweak observables beyond the one-photon-exchange approximation, one has to 
include more complicated objects. 

A very simple consistent definition of a(/x) that avoids these drawbacks is to take 
Eq. ([3]) at face value and to install Gp as an input parameter via Eqs. ([HD and (fTUli . viz. 


a(/x) 


V2 GpM 2 , 

7T 


[1 + Swill)] 


Mw 1 + Swjli) 
M 2 Z l + Szili) 


(31) 


where Swifi) and Szifi) are given in Eqs. (JHD and ([TO . respectively. This definition of 
a in) is exact and gauge independent. Since Swin) an d Szin) themselves depend on a in), 
Eq. (l3ip provides an implicit definition of a in), which may be solved iteratively, e.g., using 

7 In this work, we actually only need a(/z) through order 0{a). 
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the Newton-Raphson method. Working through two loops, we thus obtain numerically 
at n = M z 

a~ 1 {M z ) = 132.233 — 4/741+^512+^203, = 128.208, (32) 

OL OLOL s Op- 

where the first number on the right-hand side is the tree-level value a F = \/2G F My V {l — 
Myy/Mz)/n = 1/132.233 ... and the 0(a), 0(aa s ), and 0(a 2 ) contributions are specihed 
separately. Alternatively we can expand Eq. (13T|) as a power series and truncate it beyond 
0(a 2 ), as 

Oi{fx) = a F (l + a F X w (n) + a F a s (n)X n (n) + a 2 F X 20 (^) + ...), (33) 

with analytic coefficients X %] (ji) . Numerically, this gives 

a~\M z ) = 132.233 -464^+0491+^09^= 128.166. (34) 

OL OLOLs Ol 2 

The difference between the values in Eqs. (132]) and (1341) may be interpreted as the theo¬ 
retical uncertainty due to higher-order effects. 


3 Results and discussion 

We are now in a position to present our numerical analysis. For this, we adopt the 
following values of the input parameters from the Particle Data Group |35j : 

M w = 80.385(15) GeV, M z = 91.1876(21) GeV, M H = 125.7(4) GeV, 

M t = 173.21(51) (71) GeV, M b = 4.9 GeV, 

G f = 1.1663787(6) x 10~ 5 GeV' 2 , 

cr\M z ) = 127.940(14), a s (M z ) = 0.1185(6). (35) 

We neglect the masses Mf of the light fermions / ^ t,b, since their effects are negligible 
and do not play any role in our considerations. The mass of the bottom quark is taken 
into account in all one-loop expressions, but is neglected in higher-order terms. The 
order 0(aa s M 2 /M^) of the discarded terms is far beyond the accuracy of the two-loop 
approximations. In most cases, we consider the threshold corrections in Eqs. (I9l) - (fl2l) at 
a relatively high matching scale, of the order of M z or M t . However, in the case of the 
bottom quark, it is more natural to choose a somewhat lower matching scale, of the order of 
M b . Therefore, we consider y b and m b separately. Evolving the value of a(M z ) in Eq. (1331) 
with the help of Eq. (129|) at order O(aaj) to the scale + = M t , we obtain a(M t ) = 127.540, 
as already quoted in Section 12.41 Applying four-loop evolution and three-loop matching 
[37J to the value of af\M z ) in Eq. (]35l) . we obtain = 0.1081. 


3.1 Threshold corrections to the couplings and masses 


We first present the threshold corrections S x with x = W, Z, H, t , which we parametrize 
as follows 


1 + <5x(/d 


Q (h) yl,0 , a (h) ^(/d yl,l , ( a ( + ) \ y2,0 , 

4n x ' An An x \ An J x 


(36) 
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The coefficients l^ 1,0 , and Y 2 ' 0 may be evaluated numerically using a computer 

program to be published in a forthcoming paper [38j. In Appendix |Aj we list all the 
0{a) and 0{aa s ) coefficients, Y*’° and in their full analytic forms and the 0{a 2 ) 
coefficients, Y r 2,0 , in the gaugeless-limit approximation. For the reader’s convenience, we 
present linear interpolation formulae for the two-loop coefficients, Y^’ 1 and Y 2,0 , at the 
two most important matching scales, /i = Mz, M t . It is understood that the variables M t 


and Mh are 

to be taken in units of GeV. At /r = Mz , we have 



Y 1 ' 1 — 

-0.277897 {M t - 173.2) 

+ 0 {M h - 125.7) 

+ 27.8149, 


y'2,0 _ 

I W ~ 

4.83906 {M t - 173.2) 

+ 1.66627 (M h - 125.7) 

+ 2740.81, 


Y 1 ’ 1 - 
1 z ~ 

-2.15581 {M t - 173.2) 

+ 0 {M h ~ 125.7) 

- 180.9, 


V 2 ’ 0 - 

1 z ~ 

-35.6751 {M t - 173.2) 

+ 1.12839 {M h ~ 125.7) 

+ 729.793, 


ytl _ 

1 H — 

-85.325 {M t - 173.2) 

+ 67.8568 {M h - 125.7) 

- 3770.21, 


W 2 ’ 0 - 
1 H — 

-2487.1 {M t - 173.2) 

+ 917.155 {M h - 125.7) 

- 54051.1, 


= 

-5.09107 {M t - 173.2) 

+ 0.800608 {M h - 125.7) 

- 166.005, 


y'2,0 _ 

98.6532 {M t - 173.2) 

+ 5.47244 {M h - 125.7) 

+ 3318.95. 

(37) 

At /i = M t , 

we have 




*w ~ 

-0.277897 {M t - 173.2) 

+ 0 {M h - 125.7) 

+ 96.8963, 


Y 2 ’ 0 — 

18.8904 {M t - 173.2) 

+ 1.2212 {M h - 125.7) 

+ 4714.29, 


y i- 1 - 
1 z ~ 

-2.15581 {M t - 173.2) 

+ 0 {M h - 125.7) 

- 121.818, 


V 2 ’ 0 - 

1 z ~ 

-22.9901 {M t - 173.2) 

+ 0.67285 {M h - 125.7) 

+ 2450.17, 


y 1 -1 - 

1 H ~ 

-93.1689 {M t - 173.2) 

+ 40.2044 {M h - 125.7) 

- 1933.13, 


Y 2 ' 0 — 

1 H ~ 

-407.022 {M t - 173.2) 

+ 28.456 {M h - 125.7) 

- 6753.38, 


yU = 

-0.304831 {M t - 173.2) 

+ 1.02857 {M h - 125.7) 

- 79.4078, 


Y 2 ' 0 — 

I t ~ 

19.6085 {M t - 173.2) 

+ 1.64258 {M h - 125.7) 

+ 699.287. 

(38) 

We now present the MS to pole 
Specifically, we write 

mass relationships, in forms similar to Eq. 

m- 


m|(/i) , cx{fj,) i 0 

7.^-0 — -L i . A o + 

M% 4tt b 

<*(/*) a s{n) yl.l f «0)V 

47T 47T B \ 47T / 

x5° + --- 

(39) 


for the bosons B = W, Z, H and 


M f 


1 + 


a (+) ^ 1,0 , a (+) a s{v) v l,l 


47T 


-X ’ u + 


47T 47T f 


X 1 / 1 + 



2 

■y 2,0 

x f 


+ ... 


(40) 


for the fermion / = t. The coefficients X*’°, X*’ 1 , and X 2,0 may be evaluated numerically 
using the computer program to be released [38]. Convenient interpolation formulae for 
the two-loop coefficients, X *’ 1 and X 2,0 , at /i = Mz read 


X# = 43.129 {M t - 173.2) - 37.6735 (M H - 125.7) + 2287.33 , 
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■y'2,0 _ 

v = 

\^ 2,0 _ 

A z — 

yl,l _ 

A H ~ 

A H ~ 

X}' 1 = 

X 2 ’ 0 = 


(M t - 

- 173.2) 

- 727.366 

(M h 

- 125.7) 

+ 

33860.5, 

(M t - 

- 173.2) 

- 37.6735 

{M h 

- 125.7) 

+ 

2078.62, 

(M t - 

- 173.2) 

- 704.11 

{M h 

- 125.7) 

+ 

31005.9, 

{M t - 

- 173.2) 

+ 30.1834 

{M h 

- 125.7) 


- 1510.7, 

(M t - 

- 173.2) 

+ 636.554 

{M h 

- 125.7) 

— 

31764.7, 

(M t - 

- 173.2) 

- 18.1532 

{M h 

- 125.7) 

+ 

967.983, 

(M t - 

- 173.2) 

- 352.559 

{M h 

- 125.7) 

+ 

18890.7. 

and 

flUD, we 

observe that 

the coefficients X 

■i,i 

, t and X } 


( 41 ) 


w,z,t 


are 


much larger than their counterparts Y^ zt and Y V[ } zt , typically by an order of magnitude. 
This is due to the tadpole contributions that dominate the former coefficients, but largely 
cancel in the latter, as already discussed in Subsection 12.21 In the case of the Higgs boson, 
however, the situation is quite different. In fact, Xjyand X^’ 0 ^ are in the same ballpark 
as y£’ 0) and Y^ fi \ respectively, because the tadpole cancellation does not take place. 


Table 1: QCD, 0{a), 0{aa s ), and 0{a 2 ) contributions to m t (M t ) — M t in GeV for M H = 
124,125,126 GeV. The QCD contribution includes the orders 0{a n ) with n — 1,2, 3. The 
numbers in parentheses are obtained in the gaugeless-limit approximation. 


M h [GeV 

QCD 

0 («) 

0{aa s ) 

0{a 2 ) 

total 

124 

-10.38 

12.08 

-0.39 

-0.99 (-0.47) 

0.32 

125 

-10.38 

11.88 

-0.39 

-0.96 (-0.45) 

0.14 

126 

-10.38 

11.67 

-0.38 

-0.94 (-0.44) 

-0.03 


Let us now consider the various contributions to the mt{M t ) — M t shift. In Ref. [IT], 
we presented a detailed analysis thereof, providing the contribution of order 0{a 2 ) in 
the gaugeless-limit approximation and comparing it with the well-known contributions of 
orders (D(a ™) with n — 1, 2, 3, 0{a), and 0{aa s ). Our results were summarized in Table 1 
of Ref. [IT]. Here, we improve this analysis by including the full 0{a 2 ) contribution and 
updating the input parameters as specified in Eq. (|35lh Our new results are presented 
in Table [D We observe that the gaugeless-limit approximation successfully predicts the 
true sign of the 0{a 2 ) contribution, but significantly underestimates its magnitude by 
accounting for only about one half of it. The seemingly poor quality of this approximation 
may be partly ascribed to the presence of a heavy particle on the external lines of the 
two-point diagrams involved. On the other hand, the absolute deviation by itself may be 
considered acceptable. 


Table 2: Same as in Table [1] but for S t {M t ) in units of 10 4 . 


M h [GeV] 

QCD 

0(a) 

0{aa s ) 

O(X) 

total 

124 

-599.3 

13.5 

-4.4 

2.7 (3.1) 

-587.4 

125 

-599.3 

13.2 

-4.3 

2.7 (3.1) 

-587.7 

126 

-599.3 

12.9 

-4.2 

2.7 (3.1) 

-587.9 


In Table [2] the analysis of Table [Tjis repeated for 5 t {M t ). As expected by reason of 
the tadpole cancellation, the QCD corrections are by far dominant. In fact, the elec- 
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troweak corrections only make up approximately 2% of Furthermore, we observe 

that the electroweak perturbative expansion exhibits a useful convergence behavior. In 
contrast to the case of the ) — M t shift, the gaugeless-limit approximation works 

well, overestimating the true 0(a 2 ) correction by about 15%. 

Table 3: Same as in Table [I] but for 5n{M t ) in units of 10~ 4 . There are no pure QCD 
corrections in this case. 


M h [GeV] 

OH 

0(aa s ) 

OH) 

total 

124 

-114.8 

-107.5 

-26.6 (-29.1) 

-248.7 

125 

-114.5 

-105.2 

-26.4 (-29.2) 

-246.1 

126 

-114.1 

-103.1 

-26.3 (-29.3) 

-243.5 


In Table [3[ the corresponding results are presented for A (M t ), Of course, there are 
no pure QCD contributions in this case. As already observed in Ref. [16], the 0(aa s ) 
contribution almost doubles the 0(a ) one. The 0(a 2 ) contribution reaches about one 
quarter of the (D(a ) one. The gaugeless-limit approximation works here even better than 
in the case of S t (M t ), with a deviation of about 10%. 

Finally, we study the quantities related to the bottom quark, at n = M b . For the MS 
to pole mass relationship, we obtain 

{m b (M b ) - M fe } QCDi0(Q)i0(aQs)i0(Q2) = -0.85 - 1.90 - 1.53 + 1.75 (1.80)GeV. (42) 

In Eq. (l42lh the electroweak corrections are overwhelming and do not exhibit a useful 
convergence behavior, again because of the uncanceled tadpole contributions. The latter 
render the electroweak extension of the MS mass definition for the bottom quark quite 
unfeasible in practice, while they do not affect its pole mass. This situation is unfamiliar 
from pure QCD, which is tadpole free. Here, the MS mass is frequently preferred to the 
pole mass because it is only sensitive to short-distance effects, while the latter suffers 
from a renormalon ambiguity. This is reflected by the relatively large size of the Erst 
term on the right-hand side of Eq. (142|1 and the slow convergence of the 0(a ”) corrections 
with n — 1,2,3 which build it up. We also learn from Eq. (142|) that the gaugeless-limit 
approximation works here remarkably well, within an error of about 3%. On the other 
hand, the threshold corrections to the Yukawa coupling, 

{1 + Sb(M b )}Q C D,o(a),o(aa s ),o(a 2 ) = 1 - 0.1728 - 0.0190 - 0.0112 + 0.0032(0.0033), (43) 

are perturbatively stable as for the tadpole contributions. Also here, the gaugeless-limit 
approximation works at the 3% level. 

4 Conclusions 

In this paper, we analytically evaluated, through two loops in the SM, the threshold 
corrections to the electroweak gauge couplings, the top and bottom Yukawa couplings, 
the quartic self-coupling, and the vacuum expectation value of the scalar field as well as 
the MS to pole mass shifts of the top and bottom quarks. Specifically, we included the 
corrections of orders O(a), 0(aa s ), and 0(a 2 ). We emphasized the importance of the 
tadpole contributions to render these corrections gauge independent. Besides comparisons 
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with the literature, UV finiteness and RG invariance served as further checks of the 
correctness of our results. 

The threshold corrections to the gauge and Yukawa couplings are finite in the limit 
of vanishing Higgs-boson mass Mp due to cancellations of leading tadople contributions, 
and their perturbative expansions exhibit useful convergence. The threshold correction 
hf/(/r) to the quartic scalar self-coupling A(/r) scales as 1/M)/ for M H —> 0, but A(/i) is 
finite in this limit. Also in this case, perturbative stability is intact. All these threshold 
corrections are central ingredients for RG analyses within the SM and, in particular, for 
the determination of the Mp lower bound from the requirement of vacuum stability at the 
scale of the Planck mass [I 6 J UHl ED] • By contrast, the MS to pole mass shifts, which do not 
enter such RG analyses, suffer from sizable tadpole contributions, which are particularly 
severe in the case of the bottom quark. As a way out of this problem, it was proposed in 
Refs. ITS], ITT] to define the running fermion masses in terms of the MS Yukawa couplings, 
as m Y j(n) = 2~ 3,/4 Gp 1/2 y f (n). 
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A Analytical results 

In this appendix, we present analytical results for the threshold corrections A r and S x 
with x = W, Z, H,t,b defined in Eqs. (t8|) (TT2|) . The MS to pole mass relationships for 
these bosons and fermions then follow according to Eq. (TT3j) . 

Let us first introduce the notations. We work in dimensional regularization with 
d = 4 — 2e space-time dimensions and ’t Hooft mass (i. N c = 3 is the number of quark 
colors, Cp = [N'l — 1)/(2Y C ) = 4/3, and uq = 3 is the number of fermion generations. 
We express our results in terms of the pole masses and use the abbreviations 

w = M 2 W , z = M|, h = M 2 Hl t — M 2 , b = M 2 , 

Si, = 1 — w/z, L x = ln(Y//i 2 ) (x — w, z,h,t,b). (44) 


As is well known, any one-loop two-point Feynman integral can be expressed in terms 
of two types of master integrals, namely 


A(ui) 


B(q 2 -,u lj u 2 ) 




dk 

u i][( k - q ) 2 


u 2 \ ’ 


(45) 
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where 


dk = (e 7 V) £ ^ (46) 

is the integral measure in d-dimensional Minkowski space-time, Uj are the squares of the 
masses of the internal propagators, and q is the external four-momentum. To suppress 
the appearance of Euler’s constant 7 e in the Laurent expansions in e, we pull out the 
factor e~ lE£ for each loop momentum. 

In the following, we shall need the expansions of the Feynman integrals in Eq. (IT5l) in 
e through order 0(e). They read 


where 


A (tti) 


B(q 2 ;u 1 ,u 2 ) 


i — - A q (ui) - eA 0 , £ (ui) + 0(e 2 ) , 

L £ J 

i ^ + B 0 (q 2 -,u 1 ,u 2 ) + eB 0t£ (q 2 -,u 1 ,u 2 ) + 0(e 2 ) , 


Ao(wi) 

A 0 ,s( u 1 ) 

B 0 (s;ui,u 2 ) 

B 0t£ (s;ui,u 2 ) 


^(hr^-fj , 


1 

— / dx In 

0 


XUi + (1 — x)u 2 — x(l — x)s 
k 2 


1 


1. 1 f , . 9 XU 1 + (1 - x)u 2 - x(l - x)s 

2 Ci + 2.l dx hr -y- 


(47) 

(48) 

(49) 

(50) 

(51) 


The sign conventions in Eq. (l47l) have been chosen in accordance with the program library 
TSIL [39j. Some care must be exercised when a complex-valued function B 0 gets squared 
or multiplied by some other complex-valued Bo function. In particular, one needs to 
distinguish [Rei?o(- • • )] 2 from ReRg(- • • )• In such cases, we explicitly indicate when the 
real part should be taken. In particular, we have 


Bo(t] 0, 0) = Re B 0 (t ; 0, 0) — m, 


(52) 


if t > 0 , and 

t — VO 

B 0 (t] 0, w) = Re B 0 (t] 0, w) — in -, (53) 

if t > w > 0. The former situation occurs, e.g., for 5h and 5z and the latter for 5 t . 

At two loops, the most general form of the self-energy integral is given by 

Jaia2a,3a4a,5 (Q ■ 1 R, ^ 2 , ^3, ^4? ^ 5 ) 

f _ dh dk 2 _ 

J [k\ - UxYOkl - u 2 } a 0 (ki - q) 2 - u 3 ] a 3 [(k 2 - q) 2 - u 4 ] a4 [(^i ~ ~ u 5 ] a5 ’ 

where aj indicate the powers of the respective propagators. The numerator of a Feynman 
integral may be incorporated in the representation of Eq. (15T|) by allowing for some of the 
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indices dj to be negative. Alternatively, one may reduce any tensor integral to scalar ones 
in higher (shifted) space-time dimensions [40j. Then, the recurrence relations of Ref. [41] 
may be used to shift the indices as well as the space-time dimension to the appropriate 
values. After the reduction, any two-loop two-point Feynman integral may be expressed as 
a linear combination of a finite set of master integrals with the coefficients being rational 
functions of q 2 , Uj , and d. However, the choice of the master integrals is not unique. In 
particular, our set of master integrals differs from the one in Ref. [41] . A drawback of 
the set in Ref. [41 J is that the coefficients in front of master integrals may contain poles 
in e. When this happens, then deeper e expansions of the master integrals are required]))] 
Furthermore, the master integrals themselves may posses infrared divergences besides the 
ultraviolet ones. It is possible, however, to choose the master integrals in such a way that: 

(a) the coefficients in front of the master integrals always have smooth limits e —> 0 and 

(b) the master integrals themselves are infrared finite. Thanks to property (a), we do not 
need to expand master integrals beyond the order O(e 0 ). 

We keep our notations for the master integrals very close to those of Ref. [39]. Specif¬ 
ically, the different integrals in Eq. (154]) are denoted by the bold-faced letters I. S, T, 
U, V, and M. They always have indices dj = 0,1, except for T and V, which have 
dj = 2 one time. These functions generally have ultraviolet poles in e. We denote the 
O(e 0 ) terms of these master integrals by Jo, So, To, To, Vo, and M 0 , respectively]^] The 
arguments of all these functions are written as in Ref. [39], except that we always include 
q 2 as the first argument if it is present. As in Ref. [39], we exclude /i 2 from the argument 
lists. 

Specifically, J 0 , So, T 0 , To, Vo, and M 0 are defined via Eq. (154|) as 


<Tmi(g 2 ; u i, m 2 , U3, m 4 , u$) — M 0 (g 2 ; u±, u 2 , u 3 , 14, u 3 ) + 0(e ), 


JoilOliq ;Ui,U 2 ,U 3 ,U 4 ,U 5 ) — 


m 2 T M3 T W5 

2 F 2 


+ ( A 0 (m 2 ) + A 0 (u 3 ) + A 0 (u 5 ) - U2 + U 3 + u5 + O I 


+ So(cT; m 2 , M3, u 5 ) + 0(e ), 
Jo210l(q~] Ml, M2, M3, M4, M5) = + 


^0(^2) . 1 


m 2 + 2 J e 


To(q 2 ] m 2 , m 3 , m 5 ) + 0(e ), 


Jnooi (g 2 ; Mi, m 2 , m 3 , m 4 , m 5 ) — 


Mi + M 2 + M5 

2F 2 


Ao( Ml) + A 0 (m 2 ) + Aq(m 5 ) - — + 112 + U ° ^ i + J 0 (mi, m 2 , m 5 ) + O(e ), 


Jiuoi(q ;mi,m 2 ,m 3 ,m 4 ,m 5 ) — 


Mi + M 2 + M5 


2e 2 

~ ( - + B 0 (q 2 ] Ml, m 3 ) )- u 0 (q 2 ] m 3 , mi, m 5 , m 2 ) + 0(e) , 


2e 2 


8 See, e.g., the explicit calculation in Ref. [5], based on the setup of Ref. m, where poles in e through 
the second order arose and some the two-loop master integrals needed to be expanded through 0(e 2 ). 

9 The difference between these functions and their counterparts I , S, T, U, V, and M defined in 
Ref. [39] may be gleaned from Eqs. (2.34)—(2.39) therein. 
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4noi(<? 2 ; Ui, U 2 , u 3 , u 4 , u 5 ) = i(q 2 + u 1 - u 3 ) ( B 0 (q 2 ; ui, u 3 ) - 1 ) 


2Mui) + -—^oM 

U 3 

V 0 (q 2 ] u 3 , ui, u 5 , u 2 ) + 0(e ), 


A(g 2 ,Mi,M 3 ) £ 


( 55 ) 


where A (a;, y, z) = x 2 +y 2 + z 2 — 2xy — 2yz — 2zx. As in the one-loop case, the above master 
integrals may develop imaginary parts, while our threshold corrections take strictly real 
values. Since the two-loop functions never appear in products with complex coefficients, 
we impose the rule that their imaginary parts are to be discarded. 

In Ref. [IT], the threshold corrections Ar, 5 t , and 5 b were calculated at 0(a 2 ) in the 
gaugeless-limit approximation and expressed in terms of the functions B(x) and <l>(x) 
defined, respectively, in Eqs. (41) and (42) therein. The latter are related to the function 
Jo introduced above with specific arguments, namely 


Io(h,t,t ) = i(4f-/r)<f> (J^j - i(7 + C 2 )( 2 t +h)+ 6 tln-^+ 3/rln-^ 


71 ^ 1 . .. O 1 / . 7 \ - O t 

h ln — ln — - 2 h ln 75 “ 2 _ ln “ ’ 


4(0 ,M) = 


y, yz 

— n A 

2 V t 




~ (7 + C 2 ) (t + h) + 3 1 ln — + 3 h ln — 
2 y 2 y 2 


i + ^ln4 + i n l 






— h In 2 —z — t In 2 -^- 
y 2 y 2 


(56) 


Similarly, the functions R and I 2 dehned in Eq. (18) of Ref. [IT] are related to the function 
B 0 dehned in Eq. (j5U]l above, as 


= -ln+-/i (A, 

= -Ink_/ 2 ^ 1 V (57) 


In the remainder of this appendix, we list our analytic results for the threshold cor¬ 
rections Ar and 5 X with x — W, Z, H, t, b through two electroweak and three QCD loops. 
The pure QCD corrections, of orders 0(a ”) with n — 1,2, 3, only arise for and 5ft. In 
the case of 5 t , the bottom quark is treated as massless. In the case of 5b, the top quark 
is decoupled. The O(a) and 0(aa s ) corrections are exact, except that the light-fermion 
masses are neglected. In the 0(aa s ) corrections, also M b = 0 is put, except in A 0 (5). Our 
exact formulae for the 0(a 2 ) corrections are too lengthy to be presented here, so that 
we resort to the gaugeless-limit approximation. The corresponding results for 5 W and 5 Z 
vanish, which provides a welcome check for our calculation. We are thus left with those 
for A r, 5 h , 5 t , and 5 b . Our master formula reads 


1 + 5 


1 + 5qcd + 


167T 2 


X 


1,0 


+ 


9 9t 


167T 2 167T 2 


CpX 1,1 + 



(58) 


where 5 = A r,5 x and the coefficients x 1,0 , x 1,1 , and x 2,0 correspondingly carry the sub¬ 
scripts A r and x — W, Z, H, t, b. 
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A.l <5qcd 

The pure QCD correction in Eq. 


is given by 


<5qcd(a0 — 


(X. s 

Air 




+ n h 


143 


/J.M 

32( 2 


+ 


47T/ 


,71 

+ "' | y+ 


3305 

18~~ 

I 6 C 2 


64C2 + 8C3_32C2 ln2 
3 3 3 


+ 1 


/iM 


314 . 52 

■— + {n h + n l )— 


+^m ( -14 + (n h + ni)- 




1259285 99980(2 584( 3 


4864a 4 6820(4 4928( 2 , 

H-r-r— In 2 


27 

13640(5 , 10648(2(3 , 896(2 , 2 


9 


162 

+ 


608 


+ Yf ln 2 + rih 

81920(2 


9 3 27 

7315526 215728( 2 6008( 3 


V 243 


81 


27 


9 

512a 4 

27 


+ 


9 

6560( 4 

27 


9 

hT 2 


27 


ln 2 - — In 4 2 - 80( 5 + 96( 2 ( 3 + In 2 2 ) + n, ( 172318 
81 S S S 27 ) V 243 


15056(2 5656( 3 512a 4 4880( 4 1408( 2 , „ 64, , 

--—-—-“j - -In 2 — —111 2 

27 27 27 27 27 81 


In 2 2 + n 2 


256( 2 
27 

416( 2 224( 3 \ 


18962 , 512(2 , 352^ , 2 

+ + —— I + n t 


729 


135 


27 ) 


4706 

~729~ 


+ n h 7ii 


27 27 ) 

-192( 2 ln 2 - 384( 2 + 48( 3 + n h ( 4 


23668 , 416(2 , 128(3^ , , 

H-^-1-I + (m 


729 27 

18052 


42650 


714164 1120(9 448( 3 128( 2 

+ ni I _ 4--— 4---h 


4 -n, 


V 27 


712 


9 

64( 2 


-( n h + ni ) 


81 9 

104 


9 

4- n h ni 


9 

2288 

’~&T 


27 9 

In 2^ 4- n\ 
64( 2 


27 ) 

1472( 2 + 448G + 


9 

128(2 


9 

1576 


+ 


9 

128(2 \ 


ln 2 


+ 


9 


-1- l 2 

W f /zM 


81 1 9 J 

3034 


4- (n h + ni ) 


27 


/3 


, ,128 . . 2 16 
-81 4- (n h 4- n t ) —- (n h 4- n *) 2 — 


428 

IT 

(59) 


where rp is the number of massless quarks, rih = 1 is the number of heavy quarks with 
pole mass M, = ln(/i 2 /A7 2 ), and a 4 = Li 4 (l/2) « 0.517479. 


A.2 A r 

The coefficients x 1,0 , x 1 ’ 1 , and x 2,0 in Eq. 




for Ar read 
1 1 


_ _ 16 

N c +- N c ) 4- A 0 (b)( - ; - 

c 2 w ’ V A 2w(t-b) 


N c - 2-— N c ) 
hw 


,,,,,3 1 3 1 , . . , .3 a? 3 1 3 1, 

+ Mh)( J- - 4 + a ^ 2M, + llsl~ ii 1 


, , ,, 3 1 
4- A 0 (m7)(—-- 


4zS 2 


3 1 3 13 1 . 

h 2 w 4 (h — w)' 


1 z 2 1 b 2 AT 

-4- N r 

4 hw 4u>(f — 6 ) 


tb 


1 f , T 1 h 1 z w 

-N -\ _ N -4 2_ 

4w(t — b) c 4 w c 8 w 8 w h 


(60) 
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( 61 ) 


1.1 

Nc 

X Ar = 

w 

2,0 

1 

Ar,gl 

w 2 


£ 


5 

2 J 


3 h 1 11 £ 

327 + 16 + ~%h, 


t 2 37 

'J-Y 1 


5 5 t \ 

16 + 8 h) 
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Equations fl60]) -(l62l) agree with Eqs. (37)-(43) in Ref. [17] . 
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Equation (1711) agrees with Eq. (A.34) in Ref. [16] . An expansion of Eq. (T72T) in A H = 
1 — Mjj/M 2 may be found in Eqs. (19)-(23) of Ref. [17]. 
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Equations (I73l)-(l75l) agree with Eqs. (16)-(18) of Ref. [17] . Notice that Eq. (173)1 contains 
terms that behave as 1/5 for Mi —> 0. However, these singularities cancel rendering xl’° 
finite for Mi = 0. This may be seen by expanding the functions B 0 (b] z; b ) and -Bo(5; w, t ) 
through order 0(b), which yields 
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B Higgs-boson mass renormalization constant in the 
MS scheme 


The evaluations of the threshold corrections S x with x = W, Z, H, t, b at two loops, as 
described in this work, require the renormalization of the masses m x at two loops. In 
the MS scheme, the bare masses mp tX are expressed in terms of the renormalized masses 
m x (n) as 

m 0,B = ' m l(h)^(h), B = W,Z,H , 

moj = m /(h)Z/(h) ; f = t,b, (77) 


where Z x are the mass renormalization constants. The two-loop expressions for Z and 
Zz may be found in Refs. [4J [5] and those for Z t and Z b in Ref. [T7j. In the remainder of 
this appendix, we present the last missing mass renormalization constant, Zjj, in terms 
of MS-renormalized parameters. We have 
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